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Modeling an efficient Brownian heat engine 
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National Central University, Jhongli, 32054, Taiwan 

We discuss the effect of subdividing the ratchet potential on the performance of a tiny Brownian 
heat engine that is modeled as a Brownian particle hopping in a viscous medium in a sawtooth 
potential (with or without load) assisted by alternately placed hot and cold heat baths along its 
path. We show that the velocity, the efficiency and the coefficient of performance of the refrigerator 
maximize when the sawtooth potential is subdivided into series of smaller connected barrier series. 
When the engine operates quasistatically, we analytically show that the efficiency of the engine can 
not approach the Carnot efficiency and, the coefficient of performance of the refrigerator is always 
less than the Carnot refrigerator due to the irreversible heat flow via the kinetic energy. 

PACS numbers: 05.40.Jc, 05.60.-k, 05.70.-a 



I. INTRODUCTION 

Brownian heat engine functions as transducer of ther- 
mal energy into mechanical work. It rectifies thermal 
fluctuations into a unidirectional current as long as the 

[ system is out of equilibrium. In the last few decades, 
the study of such a tiny engine has got considerable at- 
tention not only for the construction of a miniaturized 
engine that can help us to utilize energy resources at 
microscopic scales [l|, 01 but also for better understand- 
ing of the nonequilibrium statistical physics [3, 0] • The 
thermodynamic property of Brownian heat engine has 

\ been explored intensively by considering different model 
systems. Brownian heat engine working due to spatially- 
variable temperature is one of the model systems which 
has been studied in the pioneeringworks by Biittiker 5] , 
Van Kampen [6], and Landauer [7|]. After the work of 
Biittiker, the ratchet model has been the subject of sev- 
eral authors [8, 9, 10, 11, 12, 13]. Recently, we consid- 

I ered an exactly solvable model of the heat engine and 
investigated the conditions under which the model works 
as a heat engine, as a refrigerator and as neither of the 
two [l3]. Not only we exposed the energetics of such an 
engine at a quasistatic limit but we also found the ther- 
modynamic properties of the engine when it operates at 
a finite time. 

In modeling of Brownian heat engine, one crucial but 
unexplored issue is the way how to maximize the veloc- 
ity, the efficiency and the coefficient of performance of 
the refrigerator. For instance, in construction of artificial 
Brownian motors, one may need to design an engine that 
accomplishes its task as fast as possible and efficiently. 
The Brownian heat engine operates autonomously Q. 
The performance of this engine (how fast and efficiently 
can it achieve its task?) relies upon the way how the 
ingredients of the system are arranged prior to the en- 
gine operation. The purpose of this theoretical work is 
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to present one possible way of improving the performance 
of the heat engine. Earlier, the mean first passage time 
(MFPT) of a Brownian particle that walks over rugged 
sawtooth potential was studied by us and in the work 
0, [l3| using super symmetric potential approach and 
using properties of random walk on networks formulated 
by Goldhirch and Gefen [H, [l^ . The sawtooth poten- 
tial was systematically subdivided into series of barriers 
without altering the barrier height, the potential width 
and the area under the barrier. The theoretical works 
revealed the existence of an optimal barrier subdivision 
that minimizes the MFPT. In this work, we use similar 
approach and study the effect of subdividing the ratchet 
potential on the performance of the engine. We show that 
the velocity, the efficiency and the coefficient of perfor- 
mance of the refrigerator tend to increase as the ratchet 
potential is subdivided into series of barriers. 

Unlike macroscopic heat engines, the Carnot efficiency 
is unattainable for Brownian heat engines when the en- 
gines work quasistatically because of the irreversible heat 
flow via the kinetic energy [1, Q . In this work, we obtain 
a simple analytic expression for the efficiency at qua- 
sistatic limit. The analytic result reveals that the effi- 
ciency of the engine never goes to the Carnot efficiency 
at quasistatic limit. Another important but unexplored 
issue is the influence of the heat flow via the kinetic en- 
ergy on the coefficient of performance of the refrigerator. 
In the present work, we analytically show that the co- 
efficient of performance of the refrigerator is always less 
than the Carnot refrigerator when the engine operates 
quasistatically. 

The paper is organized as follows: In section II, we 
present the model. In section III, we study the depen- 
dence of the efficiency and the velocity on the model pa- 
rameters in the absence of external force. We show that 
the velocity and the efficiency attain optimum values at a 
particular value of barrier subdivision N. At qausistatic 
limit, the efficiency never goes to Carnot efficiency for 
any N as the heat transfer via the kinetic energy is irre- 
versible. In section IV, we consider the model in the pres- 
ence of external load. We find that the velocity, the effi- 
ciency and the coefficient of performance of the refrigera- 
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FIG. 1: Periodic sawtooth potential without load when the 
number of barrier subdivisions A'^ = 1. The sawtooth poten- 
tial is coupled with the hot Th and the cold Tc reservoirs. 



tor attain maximum values when the sawtooth potential 
is subdivided into series of smaller connected barrier se- 
ries. We also show that the efficiency of the engine never 
approaches the Carnot efficiency and, the coefficient of 
performance of the refrigerator is always less than the 
Carnot refrigerator at quasistatic limit. Section V deals 
with summary and conclusion. 



II. THE MODEL 

Consider a Brownian particle which walks in a viscous 
medium in a periodic sawtooth potential whose potential 
profile (see Fig. 1) is described by 



U{x) 



if -io < a; < 0; 
if < a; < Lo; 



(1) 



where [/q and Lq denote the barrier height and the 
width of the sawtooth potential, respectively. The vis- 
cous medium is alternatively in contact with the hot 
Th and the cold Tc reservoirs along the reaction coor- 
dinate as shown in Fig. 1. Using the same theoretical 
frame work [l^, the sawtooth potential is subdivided 
into series of smaller connected barrier series. For ex- 
ample, Fig. 2 shows the left and the right sides of the 
sawtooth potential shown in Fig. 1 are subdivided into 
three small steps N — "i. For single barrier step be- 
tween X2 and Xi, for simplicity, we choose Ui = 2U2 
and a = 2b where X4 ~ X2 = a + b. In general, for N 
equally spaced intervals, from the top of the barrier to 
either side, Lq and Uq are given by Lq — Na + {N — 1)6 
and Uq = NUi ~ {N — 1)U2- Such parameterization is 
physically reasonable as the barrier height, the barrier 
width and the area under the barrier remain approxi- 
mately constant as N is varied [l^, • 

The Brownian particle attains a directional motion 
when it is exposed to the potential coupled with spatially 
variable temperature. For such a system, the general ex- 
pression for the steady state current J for the Brownian 
particle in any perio dic p otential with or without load is 
reported in the work |14l | . The closed form expression for 
the steady state current J (please refer Appendix A, ref. 




FIG. 2: Periodic rugged potential without load for A'' = 3. 
The potential is coupled with the hot Th and the cold Tc 



reservoirs. 



[T3 |) is given by 



J 



-F 



G1G2 + HF 



(2) 



The drift velocity v of the particle is associated to the 
steady state current J and it is given hy v = 2JLq. 

The hot reservoir is the ultimate source of energy for 
the engine. When the engine works as a heat engine, the 
net flux of the particle is from hot to the cold heat baths. 
Hence when the particle moves from the hot to the cold 
heat baths, for any N , the particle takes Uq-\-v^Lq + f Li^ 
amount of energy from the hot reservoir to surmount the 
potential of magnitude C/q and to overcome the viscous 
drag force W7 as well as the external force of amount /. 
On the other hand, 1/2/cb(TK — Tc) amount of energy 
is transferred from the hot to the cold heat baths via 
kinetic energy when the particle walks from the hot 
to the cold heat baths. Hence the Brownian particle takes 

Qh = (C/o + fLo + IV Lo) + \kB[Th - Tc) (3) 

amount of heat from the hot reservoir. The heat flow to 
the cold reservoir is given by 



1. 



Qc = Uo- Loif + -/v) + MTh - Tc) 



(4) 



When the engine acts as refrigerator, the net flow of the 
particle is from the cold to the hot reservoirs. Note that, 
due to the particle recrossing between the hot and the 
cold reservoirs, heat is leaking from the hot to the cold 
reservoir of magnitude l/2kB(Th — Tc). This is in op- 
posite direction to the heat being taken out of the cold 
reservoir. Hence, this quantity contributes as negative to 
Qc- Thus, the net heat flow out of the cold reservoir is 
given by Qc^Uq- Lo{f + jv) - ^ksiTh - Tc). 

Not all motors are designed to pull loads and alter- 
native proposals for efficiency depend on the task each 
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motor performs. Some motors may have to achieve high 
velocity against a frictional drag. This basicaUy imphes 
that the objective of the motor is to move a certain dis- 
tance in a given time interval. For such motors (/ = 0), 
the useful work is the difference between Qh and Qc'. 
W = Qh — Qc = 2'-fvLo. For motors designed to pull 
loads / ^ 0, the useful work is given by W = 2/_Lo. The 
efficiency r] and the coefScient of performance (COP) of 
the refrigerator Pref of the engine is given hy r] = W/ Qh 
and Pref = Qc/W. 

The purpose of this theoretical work, for given param- 
eter values of f/o and Lo, to find the velocity, w, the eflR- 
ciency, "q, and the coefficient of performance of the refrig- 
erator, COP, for various values of barrier subdivision, 
A''. Next, the energetics of the Brownian heat engine will 
be explored as a function of model parameters both in 
the absence and in the presence of external force. 




N = 1 



N = 2 
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III. THE EFFICIENCY AND THE VELOCITY 
IN THE ABSENCE OF EXTERNAL FORCE 

In the absence of the external force / = 0, the analyt- 
ically obtained steady state current for iV = 1 is given 

by 

where 7 denotes the coefficient of friction of the Brownian 
particle. The magnitude of the coefficient of friction of 
the Brownian particle 7 depends on temperature of the 
viscous medium. As approximation, it is considered to 
be a constant. In Appendixes A and B, the expressions 
for F, Gi, G2, and H are given for = 2 and N = 3, 
respectively. For N > 4, the expressions for F, Gi, G2, 
and H are lengthy and will not be presented in this work. 

When one omits the heat exchange via kinetic energy 
(neglecting the term l/2kBiTh-Tc) in (3) and (4)), the 
efficiency goes to Carnot efficiency r]c at quasistatic limit. 
The quasistatic limit of the engine is obtained when Uq 
goes to zero. For any N, we explore the efficiency at 
quasistatic limit and find that lim 77 = (T/i — Tc)/Th, 

which is exactly equal to the efficiency of the Carnot 
engine. 

The heat flow via the kinetic energy significantly af- 
fects the efficiency of the Brownian heat engines. Next, 
considering the heat flow via the potential and the ki- 
netic energies, we explore the thermodynamic property 
of the engine. Let us introduce dimensionless param- 
eters before exploring the dependence of the velocity 
and the efficiency on different values of barrier subdi- 
vision N. We introduce scaled parameters: scaled length 
( = I, scaled barrier height u = Uo/kBTc, scaled current 
j = J/ Jo where Jq = ksTc/^L'j., scaled velocity v = 2j/^ 
and scaled temperature r = — 1. Here ks denotes 
Boltzmann's constant. For simplicity, it is considered to 
be unity. We also introduce dimensionless parameters 



FIG. 3: The drift velocity v versus the potential barrier u 
for values r=2 and £=1. In the limit w — > and w — > 00, 
the velocity w — > 0. w increases with the number of barrier 
subdivisions. The velocity v attains maximum values at a 

particular value of u. The potential u, at which the velocity 
of the particle is maximum, shifts to wards the left as the 
number of barrier subdivisions increase 
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FIG. 4: The efficiency 77 as a function of the barrier height 

u for parameter values of r=2 and 1=1. In the limit u ^ 0, 
r]rcv ^2/3 which is equal to the Carnot efficiency for the 
given parameter values. In the limit w — » 0, r)irr — » 0. When 
w — > 00, ?7 — > 0. 



Oii = rji/rji {i = 2,3, ... A^) where rn and rji are the 
efficiencies when N ~ i and = 1, respectively. 

The dependence of the steady state current or equiva- 
lently the drift velocity on the potential u can be analyzed 
by exploiting the analytically obtained steady state cur- 
rent (for instance, see Eq. (5) for the case N = 1 and, 
the expressions shown in the Appendixes A and B for the 
cases N = 2 and 3, respectively). The directional cur- 
rent is the result of spatial temperature difference along 
the ratchet potential. In the absence of the ratchet po- 
tential, the average velocity of the particle is zero, i.e.; 
the velocity v vanishes when m — > (see Fig. 3). In the 
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FIG. 5: The plot of a as a function of u for parameter values 
of T=1.5 and i=\. a is increasing function of u. For it < 29, 
a4 > Q2, a3,...Qio. This implies A'" = Nop = 4 is the optimal 
value of barrier subdivision. On other hand, when u > 29, 
A'^ = Nop = 6 is the optimal value of barrier subdivision. 



limit u — > 00, J — > as the particle encounters a diffi- 
culty of jumping the high potential barrier of the ratchet 
potential, see Fig. 3. The velocity v attains maximum 
value at a particular value of u. Note that the engine 
operates with maximum power at this particular value of 
u. The potential m, at which the velocity of the particle 
is maximum, shifts to wards the right as the number of 
barrier subdivisions increase as it can be readily seen in 
Fig. 3. Note that in the system we consider, the left 
and the right sides of the sawtooth potential are coupled 
with the hot and the cold baths, respectively. For such 
a system, positive velocity exhibits that the net flux of 
the particle is from the hot to the cold reservoirs and the 
engine operates only as a heat engine. Figure 3 shows 
that the velocity is positive for any N. 

For high potential barriers, the Brownian particle en- 
counters a difficulty of jumping the sawtooth potential 
when the background temperature is weak. Subdivid- 
ing the barrier along the reaction coordinates enables 
the Brownian particle to cross each small barrier with 
small thermal kicks and ultimately the particle crosses 
the high potential barrier within short period of time 
than the time taken by the particle when it crosses the 
smooth potential barrier. Hence subdividing the saw- 
tooth potential enhances the drift velocity v as depicted 
in Fig. 3. The possibility of enhancing the escape rate of 
a Brownian particle over sawtooth potential under spe- 
cific conditions was envisaged [igI . \IT\ . The analytical 
finding revealed that the escape rate of the Brownian 
particle is enhanced for subdivided reaction coordinate 
which qualitatively agrees with this work. 

The analytically determined efficiency 77 is plotted as 
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FIG. 6: The drift velocity v versus r for values of u—24 
and £=1. In the limit r ^ (since Th — Tc, the system is 
in thermal equilibrium), the steady state current j vanishes: 
j 0. The drift velocity i; intensifies as r and A'' increase. 

a function of u in Fig. 4 for the case = 1. When one 
considers the heat flow via the potential energy, in the 
limit w — > 0, rjrev — * 2/3 which is equal to the Carnot 
efficiency for parameter values of t=2 and £=1. On the 
other when we consider the heat flow both via the po- 
tential and the kinetic energies, rjirr ^ when u 
and u —>■ 00. This exhibits that quasistatic process may 
not be the best working condition for the Brownian heat 
engines which agrees with the claim of Hondou and Seki- 
moto 3. In addition, ry^rr attains a maximum value at 
finite value of u. The same figure depicts that rjirr < rjrev 
and riirr,rev — when u — > 00. 

The plot of a as a function of u is displayed in Fig. 5. 
The enhancement in the efficiency is high when r is small. 
For u < 29, a4 > a2, a3,...aio- This implies = Nop — 
4 is the optimal value of barrier subdivision, at which the 
efficiency (velocity) is maximum, for a given parameter 
values. Note that the optimal barrier subdivision Nop is 
sensitive to the choices of model parameters. For instance 
when u > 29, A^ — Nop = 6 is the optimal value of barrier 
subdivision. 

The net flux of the particle strictly depends on the tem- 
perature difference between the hot and the cold baths. 
In the limit r ^ (since Th = Tc, the system is in ther- 
mal equilibrium), the steady state current j vanishes, i.e.; 
j — > (see Eq. (5) for the case A^ = 1 and the expres- 
sions shown in the Appendixes A and B for the cases 
N = 2 and 3, respectively). The unidirectional current j 
is due to the non homogenous temperature proflle along 
the ratchet potential as the particle in the hot bath can 
easily crosses the potential barrier of the sawtooth po- 
tential than the particle in the cold bath. When the 
magnitude of the rescaled temperature r steps up, the 
tendency of the particle in the hot bath to reach the top 
of the ratchet potential hill increases than the particle in 
the cold reservoir. This leads to an increase in the current 
j or the drift velocity v as shown in Fig. 6. The same 
figure shows that the net flux of the particle intensifles 
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FIG. 8: The drift velocity ii versus the rescalcd load A for 
values u = 12, £ = 1 and t = 1. For A < 4, the net particle 
flow is from the hot to the cold reservoirs. The drift velocity 
V = when A = 4. Current inversion occurs when A > 4. 



FIG. 7: The plot of a as a function of r for values of m=24 
and £=1. 



with N. The plot of a as a function of r is displayed 
in Fig. 7. Significant enhancement of the efficiency is 
observed when r is small. 



IV. THE EFFICIENCY, THE VELOCITY AND 
THE PERFORMANCE OF THE REFRIGERATOR 
IN THE PRESENCE OF EXTERNAL FORCE 

In the presence of external force, the net flow of the 
particle depends on the magnitude of the external force. 
For large load, current reversal may occur and this in- 
dicates that the engine operates not only as a heat en- 
gine but also as a refrigerator. In the presence of con- 
stant external force /, similar to the previous section, 
the closed form expression for steady state is given by 
J = -F/{GiG2 + HF). For N = I, the expressions for 
F, Gi, G2, and H are given by 



F = 


e°--^ - 1, 




Gi = 










G2 = 


a 





- 1 



(6) 



On the other hand, H = A-\- B -\- C , where 



A 
B 

C 



Th \ a 



(a + e-"-I), 



abT^ 

1 ( Lq 



(l_e-«)(e''-I), 



V b 



{e'-l-h). 



(7) 



Here a = {Uq + fLo)/Th and b = {Uq - fLa)/Tc. The 
expressions for N > 2 are lengthy and will not be pre- 



sented in this work. The drift velocity v is related to the 
steady state current and it is given by v = 2JLo. Intro- 
ducing additional rescaled parameter: A = fLi/Tc, we 
study how the velocity, the efficiency and the coefficient 
of performance of the refrigerator behave as N varies. 

Figure 8 presents the plot of the velocity ?; versus 
rescaled load A for iV = 1, 2 and 3. As shown in the 
figure for A < 4, the load is not strong enough to reverse 
the direction of the net flux of the particle, i.e., the net 
flow of the particle is from the hot to the cold reservoirs 
and hence the model works as a heat engine. On other 
hand for A > 4, the current becomes negative and the 
model acts as a refrigerator. The particle velocity is zero 
at A = 4. In general, for any number of barrier subdivi- 
sions N the velocity v = Q, when the load is 



/o 



(t + 2)Lo' 



(8) 



When one omits the heat exchange via kinetic energy 
(neglecting the term \/2kB{Th - T^) in (3) and (4)), for 
any iV, in the quasistatic limit 0, the efficiency is 

equal to Carnot efficiency: 



lim rjc 



Th 



(9) 



and in the quasistatic limit v~ — » 0, the coefficient of 
performance of the refrigerator is equal to Carnot refrig- 
erator: 



lim Pref 



Th 



(10) 



We further investigate the thermodynamic property of 
the engine by including the heat exchange via the kinetic 
and the potential energies. At a quasistatic limit, for any 
N, the efficiency takes a simple form: 



f+->0 Ih 



(11) 
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FIG. 9: The efficiency 77 versus r when the engine works qua- 
sistatically. rjc denotes the Carnot efficiency while 77*^^ de- 
notes the efficiency of the irreversible heat engine when it 
operates quasistatically for fixed u — 2. 



where 



2Uo 



+ 0.54^ 



0.5T, 



(12) 



Note that 6 < 1 when Th > Tc. This exhibits that the 
efficiency of the engine never approaches the Carnot ef- 
ficiency at quasistatic Umit as there is irreversible heat 
transfer via the kinetic energy from the hot to the cold 
heat baths. The unattainability of the Carnot efficiency, 
at a quasistatic limit, was reported in the work Q by 
solving the Kramers equation at the stall force. Exploit- 
ing the analytical expressions (9) and (11), one can ex- 
plore how rj*^^ and rjc behave as a function of t as dis- 
played in Fig. 9. The figure shows that in the limit 
Th ^ Tc, "nlrr Vc while 77*,^ < r]c when Tc < Th- 
This exhibits that the heat transfer via the kinetic en- 
ergy is considerable when r steps up. 

In the limit v~ 0, the coefficient of performance of 
the refrigerator converges to 



lim Pref = 
V--.0 Ih 



where 



0.25^ 



0.25Tr 



(13) 



(14) 



When Th > Tc and within the region where the model 
works as a refrigerator, A < 1. This reveals that the co- 
efficient of performance of the refrigerator Pre/ is always 
less than the Carnot refrigerator at a quasistatic limit. 
The quasistatic behavior of the engine can be explored 
by exploiting the analytic expressions (10) and (13) as 




FIG. 



10: The coefficient of performance of the refrigerator 

c 



Pref versus r when the heat engine works quasistatically. P^-ec 
denotes the Carnot refrigerator while P^* ^ denotes Pref of the 
irreversible heat engine when it operates quasistatically for 
fixed It = 8. 
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FIG. 11: The drift velocity i; as a function of rescaled tem- 
perature r for fixed parameters u = i = \ and A = 2.0. 
For r < 0.5, the velocity v takes negative values and the en- 
gine works as a refrigerator within this region. At r = 0.5, the 
velocity u = while when r > 0.5, the net particle current is 
from the hot to the cold reservoirs. As shown in the figure, 
\v\ intensifies as the number of barrier subdivisions increase. 



< 



pG 



shown in Fig. 10. The figure depicts that P*^j 
and, in the limit Th ^ T^, P;^j P^^j. 

We next explore how v and Pref behave as a function 
of r for TV = 1, 2 and 3. The dependence of the veloc- 
ity V on the rescaled temperature r is demonstrated in 
Fig. 11 for TV =1, 2 and 3. As shown in the Fig. 11, 
when T < 0.5, the load is strong enough to reverse the 
net particle flow and the velocity is negative. On other 
hand, when r = 0.5, the temperature renormalizes the 
effect of the load and v — 0. For t > 0.5, the tem- 
perature gains strength to overcome the load and hence 
the current is positive in this region. Within the region 
where the model works as a heat engine, v strengthens 
with TV and r. In the region where the model works as 
a refrigerator, \v\ intensifies as TV increases as shown in 
Fig. 11. This is because subdividing the sawtooth po- 
tential into small barriers enables the Brownian particle 
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FIG. 12: The coefficient of performance of the refrigerator 
Pref versus r (for reversible case) for fixed u = 10, £ = 1 and 
A = 2.0. Pref is a decreasing function of t and it rises up 
with N. In the limit t — » 0.5, Pref — » 2 which is equal to the 
Carnot refrigerator for the given parameter values. 



and the coefRcient of performance of the refrigerator is 
investigated for different number of barrier subdivisions, 
N. We show that the velocity and the efficiency attain 
optimum values at a particular value of barrier subdivi- 
sion A''. In the presence of external load we find that 
the velocity, the efficiency and the coefiicient of perfor- 
mance of the refrigerator attain maximum values when 
the sawtooth potential is subdivided into series of smaller 
connected barrier series. 

Considering the heat exchange via the potential and 
the kinetic energies, we show that Carnot efficiency is 
unachievable for Brownian heat engines when the engines 
work quasistatically. Quasistatic c;onsidcration for the 
Brownian heat engines also reveals that the coefficient 
of performance of the refrigerator is always less than the 
Carnot refrigerator. 



to cross these small barriers at small thermal kicks. On 
the other hand, when r is small, the background thermal 
kick is weak for the particle to cross the smooth sawtooth 
potential barrier. 

Figure 12 depicts the plot of Pref versus r. The coefR- 
cient of performance of the refrigerator -Pre/ is a decreas- 
ing function of r and it attends a maximum value when 
N steps up. When the rescaled temperature r increases, 
Pref declines towards the Carnot refrigerator. 



V. SUMMARY AND CONCLUSION 

In this work, we consider a model of Brownian heat 
engine. The dependence of the velocity, the efficiency 



In this work, considering an exactly solvable model, 
we explore the energetics of a Brownian heat engine not 
only at quasistatic limit but also at any finite time. This 
theoretical work suggests that the performance of the 
heat engine can be improved by subdividing the sawtooth 
potential into series of small barrier steps systematically 
by considering physically reasonable parameterization. 
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APPENDIX A 

In this Appendix we will give the expressions for F, Gi, G2 and H which define the value of the steady state 
current, J^, for zero external load case when N = 2. 



F 
G2 

Ti 

T2 
T3 



Up 2Uo 

-2 - 2e=*^c -I- 2e=*^c 



Uq_ 



J{o_ l£o, 3Lo 
2e='^'> + 2e^'^h] + — ^, 
5?7o 

2!7o 



= e 

5f/o 

bUa 

3Ln"r' -2t/o I Up Up 2Uq 

-T^e -'^c + [-6Tc + 6e ^^c Tc + 3e (Tc + n)] , 
5(7o 

= Ti(r2 + r3 + r4 + T5 + T6), 



2Up 



Up 



4 Up 



£0. Up „ "^P „ iio „ "'-^0 „ 

If 3 , 1 „ (Tc + TQUp „ 2(Tc+T^)Up 



= 54e 3TcTh 



('iTc + Th)Up „ 

- 36e 3TcTh X? 



18e 



2ITc+2T^)U„ 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 
(A7) 
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36e 3TcTh 2^2^ (A8) 
T5 = -18e~55^T7i— y2_^36eTc + 3r^ 7^2^ QgSTc + T^ (7)^(7)^ (A9) 

(^C+^hJi^O 2Uo , Up 

Te = ISe '^c'^ " Th {2Tc + Th)- 18e =*-^c + (2Tc + T/,) + 

e37ta+'¥(-45T^Tc + T^) - &ThUo + 2U^). (AlO) 

APPENDIX B 

In this Appendix we will give the expressions for F, Gi, G2 and H which define the value of the steady state 
current, J, for zero external load case when N = 3. 



- u. 



I '-'0 



F = e ^1 - 1, (Bl) 

-Up 

Lne c/fl 3Ua u„ Up 3Up Up 

Gi = — ^— — (-2-2e*'^c +2e''^c +e'^c -26*^" +2e'''^'» +6^"), (B2) 
G2 = -T^TT^ie ^'^'^ ■^c>Tc + 2e^h ^-^-c Tc - 2e *'^crc) + 

2(70 

— LnT '^o 3Uo Up 

'-{Th + 2e^Th - 2e^Th - e^{Th + Tc)), (B3) 



2?7o 

H = Ti+T2{T3+Ti), (B4) 



Ti 



Lh/Tc , ^ Up_ Vp_ 3Uo_ Up_ 

" ' (-5 + 8e^"^c - 12e*^c +4e=^c +4e*^c +e^c), (B5) 



Tp'^Ti, ^'0 -'-'0 

= ^^^^ 

T3 = {-2-26^+26^+6^+6 "^^^^', (B7) 

T4 = 2e°-^^^°^TK:+'^^ + 26°-^^'^°^'^+'^' - 6e^ - 6e^ + 8e^ , (B8) 

Uo(Th+Tc) Up(,T^+2Tc) UpiST^+Tg) 

T5 = -4e -2e ^'^f'^c +4e . (B9) 

(BIO) 
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